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ABSTRACT 

Standard accretion disc model relies upon several assumptions, the most important 
of which is geometrical thinness. Whenever this condition is violated, new physical 
effects become important such as radial energy advection and mass loss from the 
disc. These effects are important, for instance, for large mass accretion rates when 
the disc approaches its local Eddington limit. In this work, we study the upper limits 
for standard accretion disc approximation and find the corrections to the standard 
model that should be considered in any model aiming on reproducing the transition to 
super-Eddington accretion regime. First, we find that for thin accretion disc, taking 
into account relativistic corrections allows to increase the local Eddington limit by 
about a factor of two due to stronger gravity in General Relativity (GR). However, 
violation of the local Eddington limit also means large disc thickness. To consider 
consequently the disc thickness effects, one should make assumptions upon the two- 
dimensional rotation law of the disc. For rotation frequency constant on cylinders 
r sin 9 = const, vertical gravity becomes stronger with height on spheres of constant 
radius. On the other hand, effects of radial flux advection increase the flux density 
in the inner parts of the disc and lower the Eddington limit. In general, the effects 
connected to disc thickness tend to increase the local Eddington limit even more. The 
efficiency of accretion is however decreased by advection effects by about a factor of 
several. 

Key words: accretion, accretion discs - relativity - X-rays: general - radiation: 
dynamics 


1 INTRODUCTION 

Stand ard disc accretion introduced by IShakura fc SunvaevI 
successfully explain s the thermal com ponent of the 
spectra of X-ray binaries (iDone et al.l 1200^ 1. It is also ap¬ 
plied to explain the “blue bumps” of the spectral of active 
galactic nuclei but h ere the success o f the model is less cer¬ 
tain (see for instance lLawrenc idioTi)). Standard disc model 
relies on several assumptions, the most important of them 
connected to the “thinness” of the disc meaning that the mo¬ 
tion of the accreted matter does not differ significantly from 
Keplerian rotation in a plane passing through the accretor. 
Equations of vertical and radial structure may be decoupled 
then due to large differences in spatial scales and gradients 
in these two distinct directions. Thinness also implies “lo¬ 
cality” condition meaning that all the energy dissipated at 
some radius is radiated away at the same radius (because 
vertical radiation transfer is much faster). 

The basic assumptions of the thin disc model are vi- 
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dated at low and high mass accretion rates. First is con¬ 
nected with the c ooling time scales t hat become very large 
for low densities (iMever et al.ll2000bll . At large mass accre¬ 
tion rates, on one hand, the cooling timescales again become 
high because of the very high optical depths, while on the 
other hand, the radiation flux generated inside the disc be¬ 
comes strong enough for the accretion disc to approach its 

Eddington limit. _ _ 

When lEddingtonl ll 192511 introduced the luminosity 
limit, he adopted spherical symmetry. But the geometry of a 
real astrophysical system may be different. The physics un¬ 
derlying the Eddington luminosity is force balance between 
gravity and radiation pressure determined primarily by the 
radiation flujIB- Above a certain photosphere belonging to 
a star or an accretion disc, the limiting value for radiation 
flux is: 

EEdd = —g, (1) 

yc 

^ There are uncertainties connected to limb darkening and radi¬ 
ation collimation that we do not consider in this paper. 
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where >c is opacity (usually set to Thomson electron scatter¬ 
ing opacity that equals h ~ 0.34cm^ for Solar metallic- 
ity), c is the speed of light and g is the gravity in the direction 
normal to the photosphere. Spherical symmetry in classical 
mechanics implies F — L/AtvR^ and g = GM/B? and hence 
there is a luminosity limit independent of distance: 


AtvGMc 38 M _i , . 

^Edd =-~ 1.4 X 10 Ti—ergs . (2) 

XT Mq 

If the luminosity is released by accretion processes, the up¬ 
per limit for luminosity implies an upper limit for the mass 
accretion rate (as radiation pressure opposes matter infall) 
and hence: 


Mmax - 

V 


1 ^Edd 


(3) 


where g is accretion efficiency (radiation energy released per 
unit accreting mass). Again, such a limit assumes spherical 
symmetry. 

On the other hand, in accretion discs, gravity (compen¬ 
sated by centrifugal force in radial direction) and radiation, 
both directed roughly vertically, depend on both the verti¬ 
cal and radial coordinates in the disc hence the radiation 
pressure limit should be applied locally. Vertical gravity be¬ 
comes stronger with height, hence increase in local energy 
release may be compensated by increased disc thickness (we 
will hereafter use half-thickness H) up to a certain thick¬ 
ness comparable to the local radial coordinate {H ~ JJ) 
where the thin-disc approximation breaks down. Above this 
accretion rate value that we will call “critical”, accretion 
disc structure is significantly distorted inside the outermost 
rad ius where the disc become s thick (spherization radius, 
see lShakura fc SunvaevI lll973l f Section IV). Radiation pres¬ 
sure becomes most important in the inner parts of the disc 
where one should take into account general relativity (GR) 
effects to estimate rigorously the real limit when standard 
disc accretion approximation is violated. In this paper, we 
attempt to estimate accurately where does a radiatively ef¬ 
ficient relativistic disc become super-critical and find that 
the Eddington limit is scaled up by a factor of several due 
to geometrical and GR effects. 

The article is organized as follows. In Section[21 we con¬ 
sider the spherization radius and estimate the critical mass 
accretion rate in the thin-disc approximation taking into ac¬ 
count relativistic corrections. In Section [3] we consider the 
different effects connected to disc thickness that appear to 
make the effective Eddington limit even higher and besides 
change some of the properties of near-standard discs below 
the Eddington limit. In Section[4l we discuss the limitations 
of our approach and make conclusions. 


2 THIN DISC APPROXIMATION 

Thin accretion disc as introduced by IShakura fc SunvaevI 
was historically the first accretion disc model suc¬ 
cessfully applied to explain the observational properties of 
the discs in X-ray binaries and active galactic nuclei (AGN). 
Fu lly relativistic thin accret ion disc model was introduced 
bv iNovikov fc Thornel (Il973l l. Vertical radiation flux in co¬ 
moving frame: 

_ 3 GMM Q _ 3 A rh Q . . 


. Rc^ . Mxc 

where r = „,, is normalized radius, m = - - = — 

GM iEdd A-kGM 

is normalized mass accretion rate, and calligraphic letters 

correspond to the coeffici ents introduced and calcula ted by 
I Novikov fc Thor^ ll 19731 ) and IPage fc Thornd i 


B = l + - 


. 3/2 ’ 


(5) 


C = 1- - + 

^ ^ ^3/2 ’ 


( 6 ) 


Q = -^= (Vr- In - Ai - A 2 - As') , (7) 

V V 4 Tin / 

where: 

A 3 ( 0:1 -af ]„ f ] (H) 

^ a;i(a;i — a;2)(o;i — 0:3) yyrhT —*1/’ 

xi, 2,3 are solutions of the equation — 3x + 2a = 0, and 
A 2 and A 3 are obtained from Ai using cyclic permutation 
1^2^ 3. 

Radiation pressure is o pposed by the vertical gravity 
that was first calculated bv iRiffert fc Heroldl lll995l ) in the 
z/R<^ 1 limit. Co-moving vertical gravity: 


where: 


_ GMz Gr 

- ^3 ^ , 


G, = l- 


4o 3a^ 

^ 3/2 3 ^3”' 


(9) 


( 10 ) 


For a thin disc, vertical gravity is always proportional to 
the vertical coordinate 2 . At some height radiation pressure 
becomes balanced by the vertical gravity: 


3 . GM Q VC 
2 ”^ c2 B CA 


( 11 ) 


This quantity has the physical meaning of the height of a 
radiation pressure supported accretion disc. It is also useful 
as an estimate for the so-called spherization radius at which 
the thin disc approximation breaks down. Since the vertical 
gravity approaches maximum at 2 ~ R, this puts also the 
local Eddington limit for the mass accretion rate. A fully 
relativistic expression for the spherization radius may be 
given in implicit form: 


3 . 

rsph = -m ■ 


Q(^sph5^) V^(rsph;^) 

I^(^sph,u) G7-(rsph,u) 


( 12 ) 


This expression may be solved numerically for rgph. The load 
of coefficients on the right-hand side approaches unity at 
large radii, but the approach is slow enough to produce 

3 

strong deviations from the classical value of rsph = -^ 1 x 1 

given by IShakura fc SunvaevI (|l973h . These coefficients in¬ 
corporate both relativistic effects and the inner boundary 
condition at r = rin (at the last stable orbit radius). It is in¬ 
structive to compare expression GH) with its non-relativistic 
analog ue taking into account only th e inner boundary con¬ 
dition llAboImasov fc Shakurair2012al ): 


In this approximation, spherization radius exists only if rh > 
rhcr = 4.5rin In Fig. [T1 we compare these two estimates. For 


3.4 2 

^sphjNR — 2^3 
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Figure 1. Spherization radii calculated according to relativistic 
(solid lines) and non-relativistic (dotted) formulae. Thick dots 
mark the critical accretion rate for the relevant Kerr parameter 
value. The curves and dots were calculated for three different Kerr 
parameter values: a = 0 (red, middle curves), a = 0.9 (assum¬ 
ing co-rotation; blue, rightmost curves) and a = —0.9 (counter¬ 
rotation; green, leftmost curves). 

accretion rates m < 1000, relativistic corrections alter the 
spherization radius by about 20%. In Fig. [TJ solid curves 
correspond to the GR case, and dashed to non-relativistic 
approximation. In both cases, middle curves correspond to 
the Schwarzschild case (a = 0) and the right and the left 
curves - to a moderate Kerr case with a = 0.9 assuming 
co-and counter-rotation, correspondingly. 

Equation (HI has no solutions if mass accretion rate 
is small enough. If mass accretion rate is large enough, 
there are two solutions, the larger of which is interpreted 
as the spherization radius. The minimal m value where the 
equation has a solution has the physical meaning of critical 
mass accretion rate. While the spherization radius changes 
slightly due to relativistic effects, the critical mass accre¬ 
tion rate defined this way increases by about a factor of 
3. In Fig. H we compare the critical mass accretion rates 
defined straightforward as LBdd/» 7 (fl)c^, and non-relativistic 
and relativistic critical mass accretion rates defined through 
local radiation pressure balance condition. Here, the dotted 
line shows the limiting mass accretion rate in the spheri¬ 
cally symmetric case, equal to rj' ^(nn , a). The dashed curve 
corresponds to a thin Newtonian disc while the solid curve 
takes into account GR effects. 

3 CORRECTIONS FOR DISC THICKNESS 
3.1 The issue of disc thickness 

As real accretion discs become thick near the Eddington 
limit, it is instructive to check whether neglecting accretion 
disc thickness introduces any important corrections into the 



Figure 2. Critical mass accretion rate as a function of Kerr pa¬ 
rameter a (negative a corresponds to counter-rotation). Dotted 
line is the critical accretion rate for spherically symmetric accre¬ 
tion (r?“^(riu, a)), dashed line is the maximal mass accretion rate 
in Newtonian approximation, and the solid line corresponds to 
the limit in full GR. 


vertical balance. Disc “thickness” is in fact a complex of sev¬ 
eral physical effects connected to violation of the basic as¬ 
sumptions of the thin-disc approximation, including: 

(i) geometrical thickness and, as a consequence, non¬ 
linear dependence of the vertical gravity on height, 

(ii) non-Keplerian rotation and differential rotation with 
vertical coordinate that may lead to additional angular mo¬ 
mentum transfer in vertical direction; this also affects dissi¬ 
pation and energy release inside the disc, 

(iii) violation of local energy release approximation due 
to both advection and radial flux diffusion, 

(iv) existence of radial and meridional velocity fields. 

Here, we will try to approximately account for the first, 
the third and the partly the second correction. We will 
assume constant rotation frequency on cylindrical surfaces 
•uj = r sin Q = const, where 9 is polar angle. This is close 
to the rotation law expected in barotropic case where ro¬ 
tation frequency should be constan t on the so-called von 
Zeipel cylinders llAbramowic 3 I 1971 I) . In Appendix [Al we 
show that the shapes of von Zeipel cylinders are indeed very 
close to cylinders in Boyer-Lindquist coordinates practically 
everywhere except inside the ergosphere. Accretion disc en¬ 
ters the ergosphere only for very large rotation parameters 
a > 0.94. Outside the ergospheric region, for a fixed cylindri¬ 
cal coordinate w = rsinS, rotation frequency always grows 
with height that implies faster rotation (and, consequently, 
stronger effective gravity) at a surface for constant spherical 
coordinate r. Hence using coordinate cylinders zu = const 
we underestimate the growth of the gravity with height. 
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Table 1. HARM model parameters, including the parameters 
of the initial Fishbone-Moncrief torus, rmax (radius of maximal 
specific enthalpy) and (inner radius) are also given. The time 
ranges (in GMjc?) used for averaging are also given. 


Model 

resolution 

a 

^max 

fin 

number 
of loops 

time range 

Blh 

192x320 

0.9 

8.4 

4.8 

1 

2000-5700 

A50s 

128x192 

0 

19.7 

12.8 

50 

10000- 

15000 

Dl 

128x192 

0.99 

10 

4 

1 

2000-5700 


Surprisingly, in all our models, not only the lines of con¬ 
stant angular velocity are close to cylinders vj = rsinS = 
const but the rotation law between r ~ 2 and r ~ 10 con¬ 
forms with Keplerian rotation with about 3—4 per cent accu¬ 
racy (relative standard deviation from the Keplerian law) for 
9 > 45°. For larger distances from the equatorial plane, ro¬ 
tation becomes strongly sub-Keplerian. On the other hand, 
in the equatorial plane, Keplerian rotation holds to about 
1 per cent. In Fig. [3] (left panel), lines of constant angular 
frequency (fl = are shown and compared to the Ke¬ 
plerian rotation law (Hk = -|-a^), and the lines of 

constant angular frequency are shown to conform well with 
the cylinders vj = const. Keplerian rotation holds even in¬ 
side the last stable orbit. 


3.2 Numerical simulations 

Numerical simulations of thick MHD discs support the 
assumption of Keplerian rotation in the inner parts of 
thick accretion dis cs and even inside the last stable orbit 
dPenna et al. I [2^. We have considered the res ults of a 
thick disc MHD simulation from I Abolmasov! ||2014|) . namely 
the run Blh presented in this work. Simulation was per¬ 
formed using the numerical code H ARM2D that solves 
ideal MHD equat ions in Kerr metric (iGammie et al.l [20031 : 
iNoble et al.ll200a ). The main goal of the simulation was to 
form a moderately thin disc {H/R ~ 0.1) around a Kerr 
(a = 0.9) black hole and consider the motion of the matter 
inside the inner edge of the disc. As the initial conditions, 
an exa ct relativistic torus solution by [Fishbone fc Moncriej 
lll976l) . section HIb, was used with small magnetic fields in¬ 
troduced that are subsequently amplified by MHD instabil¬ 
ities in the disc. Amplified magnetic fields lead to outward 
angular momentum transfer and consequently to accretion 
from the torus through a moderately thick disc formed dur¬ 
ing the first two thousand dynamical time-scales {GM/c?) 
of the si mulation. We also u se the simulation run A50s pre¬ 
sented in lAbolmasovI ll2014h that was calculated for a = 0 
and multi-loop initial magnetic field and uses smaller reso¬ 
lution. Main details of the simulations used in the present 

paper are given in Table [T] _ 

_ A thicker torus in the model of [Fishbone fc Moncriej 

([l976l) is also more extended in radial direction. As the re¬ 
sult, it may become impossible to maintain the boundary 
condition at the outer boundary, and the simulation becomes 
unstable. A way around is to consider a black hole with a 
higher Kerr parameter. More rapid rotation allows to pro¬ 
duce a thicker disc without making the initial torus model 
extended beyond t he simulation dom ain. Apart from the 
models described in [Abolmasov[ |[2014l ). we have calculated 
a similar model with a higher Kerr parameter, a = 0.99. Ab¬ 
sence of heating and cooling sources in HARM2D makes it 
difficult to simulate very thick as well as very thin discs. 
While the initial torus has H/R ~ 0.5, the accretion disc 
forms from its inner cusp that has lower geometrical thick¬ 
ness. For the Blh model, the relative thickness is H/R ~ 0.1 
as measured by Gaussian fitting, while for the model with 
a = 0.99 that we will hereafter dub as Dl, the thickness was 
around 0.2. In fact, all the equilibrium torus models, how¬ 
ever thick, h ave smaller thickne sses near the inner rim (see, 
for example, [Penna et al.[ ([2013[ ). Fig. 3). 


Our results are at odds with so me of the comprehen - 
sive thick-disc simulations such as bv [Sadowski et all ([2014f ) 
that, for a similar set-up, find rotation sub-Keplerian by 
not less than 15 per cent, both outside and inside the last 
stable orbit ( see F ig. 11 in this work). On the other hand, 
[Shafee et al.l l[2008[) . for instance, find the angular momen¬ 
tum close to Keplerian just outside the last stable orbit. 
The overall be haviour of the net a ngular momentum in the 
simulations by [Shafee et al.[ ([2008[ ) conforms well with the 
conventional picture of hydrodynamic transonic accretion: 
nearly Keplerian rotation outside the last stable orbit, Ke¬ 
plerian rotation at the sonic surface and roughly constant 
angular mo mentum (Uui] in the s upersonic region. On the 
other hand. l^dowski et al.[ |[2014h find the net angular mo¬ 
mentum rapidly decreasing inside the last stable orbit. 

The reason for this discrepancy lies probably in the di¬ 
verse structure of magnetic field fo rmed during disc accre- 
tion with frozen-in magnetic fields |[Penna et al.[[20f5 ). The 
initial conditions include small seed magnetic field subse¬ 
quently amplified by MHD instabilities. Amplified poloidal 
magnetic fields are then advected towards the black hole 
and contribute to the accumulated magnetic flux through 
the black hole horizon. This magnetic flux grows gradually 
with time that leads to secular evolution of the deviations 
from Keplerian rotation. In simulations Blh and Dl, we as¬ 
sumed the initial magnetic field potential s caling with the 
rest-m ass density as oc p. However, in [Sadowski et al.[ 
|[2014h . the magnetic field obeys a different, more compli¬ 
cated law making it stronger further from the black hole. 
This leads to a gradual increase in magnetization of the ac¬ 
creting matt er even during quasi- stationary accretion (see 
also Fig. 6 of [Sadowski et al](2014[ ). 

We conclude that the rotational structure of the inner 
disc and especially the flow inside the last stable orbit is 
strongly affected by the magnetic fields and especially by 
the magnetic field flux accumulated inside the disc. Our re¬ 
sults should work as long as the large-scale magnetic field is 
dynamically unimportant. For the time span of the Blh sim¬ 
ulation, rotation law changes very weakly after the establish¬ 
ment of quasi-stationary accretion at t ~ 2000 GM/c^, but 
simulations containing more complicated magnetic fields like 
A50s show irregular variations due to changes in magnetic 
field strength and geometry (see Fig. [S|. Angular momenta 
given in this figure are dens i ty-we ighed vertically-averaged 
values (as in [Sadowski et aP ([2014[ )). 

In our simulations, rotation frequency is everywhere in 
the thick disc close to the following “cylindrical Keplerian” 
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Figure 3. HARM2D simulation snapshot (£ = 5600GAf/c^, Blh simulation from [Abolmaso^ ll2014h h Left panel: colour-coded is 
angular velocity Q = u^ju'^ in c? jGM units, black lines are vertical {pj = const) lines corresponding to the same values of angular 
velocity as the boundaries between different colours/shades. Right panels: angular frequency slice through the equatorial plane (thick 
red points) and for all the points with polar angle cosines |/r| < l/\/2; solid black line is Keplerian law to = ((rsin g)3/2 _|_ 1 Lower 

right panel shows relative deviations from the cylindrical Keplerian law. All the distances are in GMjc^ units. Vertical dotted lines mark 
the radius of the last stable orbit. 
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Figure 4. HARM2D simulation snapshot (i = 5400 GM/c®, D1 simulation with a = 0.99). Notation is the same as in the previous 
figure. 


law: 

^{r,0)= 3,,^, . (14) 

Estimating the two-dimensional rotational velocity field is 
crucial for the vertical structure of the disc since the vertical 
structure of the disc is determined by the deviations from 
the balance between inertial and gravitational forces, the 
former dependent upon the rotation law. There is no clear 
distinction between gravity and inertia in general relativity, 


and the vertical hydrostatic balance in accretion disc may 
be alternatively viewed as a balance bet ween pressure gradi¬ 
ents in the disc and the centrifugal force llAbramowicz et ^ 

mia). 


3.3 Corrections for vertical gravity 

Let us use spherical coordinates and consider the balance 
of radiation pressure and inertial forces (including gravity) 
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6 Abolmasov & Chashkina 




Figure 5. Radial profiles of density-weighed specific angular momenta for two simulations Blh (left) and A50s (right). Black solid 
curves are averages for all the frames of the quasi-stationary phase, green dot-dashed and red dashed curves show the averages calculated 
tor the first and second halves of the period of time that we interpret as quasi-stationary accretion (2000 — 5700 GM/c® for Blh, 
10 000 - 15 000 GM/c® tor A50s). 


along the polar angle It is more convenient to use /r = cos 9 
instead of the angle itself or the vertical distance z. Relevant 
metric coefficient is 


(15) 

Then, considering a sphere of r = const, one can express 
the “vertical” free-fall acceleration as: 

=-Q.lfi = rif,,kuu'^, (16) 

where i and k may be equal t or ip, u‘*’ = flu^, and u* is re¬ 
covered from normalization condition mu'' = — 1 (we neglect 
the radial and meridional velocities here). Vertical grav¬ 
ity should be compared to the coordinate flux component 
~ P multiplied by x/c (see below Section [3.4|l . All 

the relevant non-zero Christoffel symbols are = ^gik,m 
TtMt ~ ^4 ( 11 ^) 


where: 


8 1 a 

Cir,a) = 1 - o- + — 

o T 

2 a?{7 — 4:0?) _ o' 

+ 1 - - - 12 


2 14 

3 ^5/2 


23 16 a-^ 

yyj + 

+ 17^- 


3 

8a®(l-h3a2) 16 a'* , 38 a® 

3 ^ 11/2 ^ 13/2 3 "^' 


16 a'" 


( 21 ) 


In such an approach, for most radii and Kerr param¬ 
eter values, the effective vertical frequency becomes higher 
away from the equatorial plane if the spherical radial coor¬ 
dinate is held fixed. This effect arises from the assumption 
fl = fl(rsin0) that leads to effectively super-Keplerian ro¬ 
tation near the disc surface. If one considers a more realistic 
rotation law with fl constant on von Zeipel cylinders, the 
effect becomes even stronger, except the most rapidly rotat¬ 
ing black holes (o > 0.95) where the inner parts of the disc 
exist inside the ergosphere, where the shapes of von Zeipel 
cylinders change considerably (see Appendix |A}. In our as¬ 
sumptions, ( and thus the whole next-order correction is 
always positive. We also assume that the maximal possible 
thickness of the disc is ^max = 1/^2. 


2 ar(r^ + a=) 

r tflip — i —- , fl, 


(18) 


= -4 (a^(l - M^)^^ + - a^A{l - p?)) p = 

_ A + -I- 4a‘^A + 2Ar + 

9- 

(19) 

Substituting these expressions in the expression for up 
to second order in g: 




2 


5 C(Ga) 2\ 
r\C 2 C^B ^ J’ 


( 20 ) 


3.4 Advection 

Local energy flux in an axisymmetric disc consists of diffu¬ 
sion and advection terms, the latter directed along the local 
poloidal velocity. If the advection term is small compared 
to vertical diffusion, its overall effect is to shift downstream 
the thermal energy generated in the disc. Weak advection 
in a moderately thick disc may be approximately taken into 
account by introducing a “delayed radius” assumption that 
the flux generated at radial distance vq will be then emitted 
at a smaller radius r = ro — Sr{ro). For our purposes, this ap¬ 
proach is better than assuming a fixed fraction of liberated 
energy stored in the thermal energy of the disc aa is it often 


© 2012 RAS, MNRAS 000, [THU 














Relativistic disc Eddington limit 7 



Figure 6. Flux dependence on radius for different mass accretion rates. Fluxes are divided by the mass accretion rate hence the highest 
fluxes correspond to the smallest m where advection is smaller. The thick red curve shows the standard accretion disc model, thin black 
curves correspond to rh = 1, 10, 100 and 1000 (solid, dashed, dash-dotted and dotted lines, respectively). Left and right panels show the 
cases of a = 0 and a = 0.9. 


done when considerin g discs with advection ijNaravan fc Yil 
1 19951: ISado wskil I2OIII I , since it does not include unknown 
free parameters and can account for advection dependence 
on radius. Apart from the pure flux advection, one should 
take into account that the frame where the flux is compared 
to the vertical gravity is different from the frame (or frames) 
where the energy is generated that we would further identify 
with the co-moving equatorial frame in the disc. Radial shift 
Sr is a result of interplay between vertical energy diffusion 
and radial drift. 


St — Urtdiff 5 

where Vr is the radial velocity measured in the reference 
frame co-moving with the disc, tdiff = Htv, where H = 

r/imax, and Tv = 2^^ vertical optical depth between 

the equatorial plane and the disc surface. 


Sr ~ (23) 

The main simplification is that the flux generated at the ra¬ 
dial coordinate ro is emitted in a narrow radius interval near 
r. A reasonable assumption is also to assume all the flux gen¬ 
erated near the equatorial plane and then conserved. More 
precisely, let us consider the energy flux due to radiation 
Tt = -\- u^Ft and integrate it upon a fi = Hrasx sur¬ 

face just above the surface of the disc. Conservation of the 
energy flux implies: 


Ft V~-gdSk = const, (24) 

for a given flux tube with cross-section given by the four- 
vector dSk- Setting the cross-section tangent to the 9 = 
const surface implies 

utF‘^y/—gdSf_i = const, (25) 


or 


(26) 

where quantities with ‘(d)’ index refer to the equatorial 
plane, and with ‘(s)’ index - to the disc surface, g — —p^ is 
metric determinant (hence g(d) = —ro). Since drjd) = dr^^s) 
(see above expression (1231) for Sr), area distortion is con¬ 
nected only to the difference in factors, incorporating 

relativistic effects and spherical geometry. Co-moving flux 
at the surface of the thick disc: 


Ro) = 


^jd) p 

P(„)sin6»„W 


(27) 


All the quantities with a ‘(d)’ index on the right-hand side 
should be evaluated at ro, and 9 is the half-opening angle of 
the disc. Eddington limit in polar angle direction is written 
as: 


~F(s)^gij,tt — (28) 

where is given by expression and flux F(^s) by ex¬ 

pression dig. 


3.5 Effects of disc thickness upon the local 
Eddington limit 

The effects of non-locality (including advection) are ambigu¬ 
ous: flux increases due to geometrical effects, and decreases 
due to decreasing dependence F{r) everywhere except the 
innermost parts where the disc releases more energy than a 
standard disc. In Fig. flux dependence on radius is shown 
taking into account advection and the geometrical effects de¬ 
scribed in Section lT4l Note that for very high mass accretion 
rates (m > 10^) advection is able to prevent the disc from 
surpassing the local Eddington limit. Our approximation 
does not reproduce this regime but it can definitely appear 
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when advection becomes strong and essentially non-linear. 
This advection-dominated regime may be loosely i denti¬ 
fied with the slim-disc solution (lAbramowicz et aP 1 198a : 
Isid ows i^ |201l|) . 

In Fig. [T] critical mass accretion rates calculated with 
and without corrections for disc thickness (including advec¬ 
tion) are shown. As in Fig. dotted line corresponds to 
the spherically symmetric case, and dashed line to the thin 
non-relativistic disc. Blue line (marked with a “T” in the fig¬ 
ure) shows the thin relativistic disc case considered in Sec¬ 
tion [2] However, the limiting disc thickness is taken to be 
{H/R)ina.x = (since R is interpreted as spherical rather 

than cylindrical coordinate; this is also close to the majdmal 
polar angle cosine where cylindrical rotation holds in the 
HARM2D simulation we mention above) that decreases the 
critical rate. Red (subscribed with a “G” in the figure) and 
green (marked “A”) curves show the effects of gravity and 
advection taken separately, and the solid black curve is cal¬ 
culated taking into account both effects. The critical mass 
accretion rates are shown as functions of Kerr parameter a. 
Same critical mass accretion rates are also shown as func¬ 
tions of 1 — a in logarithmical scale that allows to show in 
detail the limit of rapid black hole rotation. 

Evidently, corrections to vertical gravity increase the 
value of the Eddington limit, at least for a < 0.96. However, 
the resulting values of the critical mass accretion rate are not 
much larger than the thin-disc estimate. While an account 
of advection and pure spherical geometry effects described 
in Section [TH in fact decreases the local Eddington limit be¬ 
cause the local flux in the inner parts of the disc increases, 
stronger vertical gravity acts in the opposite direction. In¬ 
cluding both effects simultaneously increases the local flux 
limit even more than gravity alone because at larger rh, ad¬ 
vection effects decrease the local flux. The minimal possible 
spherization radius also becomes considerably larger because 
in the inner parts of the accretion disc, advection becomes 
too strong. Summarizing, we conclude that different effects 
connected to disc thickness produce corrections of differ¬ 
ent signs, but the general result holds and becomes even 
stronger. 


3.6 Zone a disc thickness 

Physical meaning of condition (I28II is vertical balance be¬ 
tween radiation and gravity, it can also be solved for to 
estimate the vertical scale of the disc. In Fig. [S] we com¬ 
pare the vertical scales of the discs for m = 30 and different 
Kerr parameters. We find that disc thickness corrections de¬ 
crease the vertical scale in the hottest and thickest parts 
of the disc while its inner rim attains finite thickness. At 
high mass accretion rates, the thickest part of the disc shifts 
to larger distances, and the surface inside the thickest part 
approaches a roughly paraboloidal or conical shape. 

For disc thicknesses, corrections for advection were cal¬ 
culated using expression (1231) for radial shift. In reality the 
radial advection correction to the pressure inside the disc 
should be smaller because the thermal energy released in¬ 
side the disc has an immediate effect upon the local pressure 
but delayed effect upon the flux coming out from the disc 
surface. Future calculations of vertical disc structure should 
take this into account together with radial diffusion of ra¬ 
diation and other effects connected to the two-dimensional 


nature of thick accretion discs. The general effect upon the 
vertical scale is smearing the abrupt decrease in disc thick¬ 
ness near the inner radius. The shape of the disc surface 
generally becomes closer to h/r = const. 


3.7 Efficiency 


Increasing mass accretion rate limit does not necessarily 
mean larger luminosity because some part of the energy dis¬ 
sipated inside the disc is advected inwards. Since m ^ 1, 
the energy is trapped in the accreting flow even inside the 
last stable orbit. Indeed, the vertical optical depth in the 
nearly free-fall region inside the last stable orbit is 

Tz = x: f {—Uik'‘)pdX, (29) 

Jo 

where A is affine parameter increasing along the photon 
trajectory and = dx^/d\. Mass conservation implies 
E = M/4-kRu^. For a vertical trajectory, dr = 0 and d6 — 0, 
and hence, if — 0: 


Tz 


xY,E = 


Em 
2ru^ ’ 


(30) 


where E = —ut is energy-at-infinity for a unit-mass particle 
in the flow. The locally measured diffusion time is tdifr ~ 
Htz/c. Radiation diffuses in vertical direction, but in radial 
direction it is advected at the local dynamical time scale 
tdyn — c(risco — rhor), where risco and rhor are the radii of 
the last stable orbit and the event horizon of the black hole, 
and the infall velocity was set to c. Relative importance of 
the two effects is given by the time ratio 

tdiff Em H . 

- — TTT-TTT ~ (31) 

fdyn 2(riSCO rhor) R 


This means that the falling flow inside the last stable or¬ 
bit starting from moderate accretion rates m ~ 1 is indeed 
optically thick to electron scattering, especially for higher 
rotation parameters when the radial free-fall time is small. 
Optical depths to different absorption processes are likely 
smaller than unity but are sensitive to clumping (all the 
opacities that scale with density squared) or to magnetic 
field strength (synchrotron absorption) and hence are more 
difficult to constrain. Below in section ITSI we estimate the 
effective optical depth to free-free absorption. 

Since the flow is optically thick one can apply the for¬ 
malism we used in Section EH The energy released in the 
disc will be shifted downstream by: 

Ar ~ utdis ipodiE. (32) 


This allows not only to estimate the overall efficiency but 
also the corrections to the spectral energy distribution that 
we will present in Section (3.81 
Accretion disc luminosity: 


Lmax — TEdd?l(n, 7ilcr)n2cr (n) — 

~ 1.4 X 10 ^® (i^) ^(“’’^cr)?Ticr(a)ergs“^ 

For a thin disc, (see, for instance, IShaoiro fc Teukolskvl 
lll986ll Section 12.7): 


ri{a, rh 



(34) 
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Figure 7. Critical mass accretion rate as a function of Kerr parameter. Left and right panels differ only in abscissa: in the left, it is 
linear in a, in the right, it is made logarithmic in 1 — a to resolve the behaviour of limiting mass accretion rates at high mass accretion 
rates. Dotted and dashed black curves are identical to those in Fig. (2] corresponding to spherical and disc non-relativistic Eddington 
limit, correspondingly. The solid blue line (also marked as “T”) corresponds to the local thin relativistic disc Eddington limit. It goes 
somewhat lower than the solid line in Fig. [^because the accretion disc thickness is limited by /i = l/\/2 rather than ^ = 1. Green curve 
(marked by “A”) takes into account advection only, and red curve (marked by “G”) only second-order corrections to vertical gravity. Solid 
black curve accounts for both effects. Vertical dotted line marks the change of sign in C, correction at the inner rim of the disc. 




Figure 8. Disc vertical scales calculated taking into account corrections for vertical gravity and advection (solid lines). Dashed lines 
are standard disc zone a thicknesses for identical parameters. Dotted straight lines correspond to ^ = r. 


In the general case, t] = 1 — E where energy E is calculated at 
the inner visible edge of the disc. Considered in the delayed 
radius formalism, E should be taken at a radius larger than 
the last stable orbit radius. For our model thick disc with 
advection, the large optical depth of the falling flow implies 

r]{a,7h) = 1 — E{rhor + Sr). (35) 

Assuming Keplerian rotation allow s to approximate the en¬ 
ergy as (cf. IPaee fc Thornd (Il974h : our estimate hear does 
not account for radial motions that should be slow as long 
as Keplerian rotation holds): 

K = (36) 


where: 

^ = + ( 37 ) 

Efficiency starts decreasing significantly from rh ~ 10 and 
then approaches a power-law asymptotic C/rh with C de¬ 
pending upon Kerr parameter. For high mass accretion 
rates, this implies a constant luminosity of CLBdd with C 
between approximately 1.3 and 2 (see Fig. H). When mass 
accretion rate reaches the critical value and the upper lay¬ 
ers of the disc become unbound, efficiency becomes smaller 
than for the standard thin disc by a factor of about a factor 
of 5 (see Fig. [iQl). Some, not very large, part of this luminos- 
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Figure 9. Efficiency dependence on mass accretion rate calcu¬ 
lated in delayed radius approximation for different Kerr parame¬ 
ters (a = —0.9, 0, 0.9 and 0.998, shown by black solid lines, from 
bottom to top). Dotted lines are thin-disc efficiencies for corre¬ 
sponding Kerr parameters. Red dashed straight lines are 1.3/m 
and 2/m. 



a 


Figure 10. Efficiencies for a thin disc (dotted), thin disc with 
advection (dashed) and thick disc with advection (solid black and 
dot-dashed green lines). The green curve also takes into account 
the radiation escaping from inside the last stable orbit. The three 
latter curves were calculated for corresponding critical mass ac¬ 
cretion rates. 


ity is emitted from the disc extention inside the last stable 
orbit. As advection enters non-linear regime, and the disc 
starts losing its matter, a logarithmical correction to this 
quantity should appear that we will discuss in more detail 
in Section 0] 


3.8 Spectral energy distribution 

As one can see from previous sections, including thickness 
effects makes the flux generally higher for intermediate mass 



Figure 11. Temperature-luminosity plot (energy of maximal 
in the spectrum was used as a measure of temperature) for 
black hole accretion discs with different mass accretion rates be¬ 
tween 0.1 and 500 spaced evenly in logarithmic space. Dimen¬ 
sionless mass accretion rates are given as black and red numbers. 
Black circles are for a thick advection-modified disc, red open cir¬ 
cles for a thin standard disc. Solid curve corresponds to the ap¬ 
proximate relation for a standard accretion disc where L oc 
Dotted curve is an estimate for the maximal effective temperature 
expected in discs with outflows. 


accretion rates and lower for m > 100. Thus, the spectral 
energy distribution first becomes hotter, then cooler with 
mass accretion rate. We have considered spectral energy dis¬ 
tribution of a multi-blackbody disc for a fixed configuration 
(Schwarzschild blac k hole viewed at 30° in clination). Public 
code GEOKERR (see iDexter &: Agoll (|2009l) for description) 
was used to compute the geodesics, and the observed in¬ 
tensity was re-calculated using the conservation of energy 
and angular momentum along the trajectories. Details of 
in tensity transformation r nay be found e. g. in section 3.3 
of I Abolmasov &; Shakur3 (|2012lj) . Black hole mass was set 
to M = IOMq. For different mass accretion rates, spec¬ 
tral energy distributions were calculated with and without 
considering thickness effects. For every spectrum, we calcu¬ 
lated integral luminosity and position of the maximum of 
vFi, (flux per logarithmic frequency interval). The resulting 
luminosity-temperature plot is given in Fig. 1111 

Because in our model the disc retains its multi- 
blackbody nature, it does not shift significantly from the 
standard disc luminosity-temperature relation. Most of the 
luminosity is emitted within the area comparable to the area 
subtended by the last stable orbit hence the luminosity is 
always proportional to the local flux in the inner parts of 
the disc and L oc T^ax ~ where Fmax is the energy 

where the spectral energy distribution has the maximum. In 
Fig. 111! solid curve is the approximate temperature in the 
disc without correction factors: 

Tmax = 5.29 j rr^/"keV (38) 

For the set of parameters we use here, Tmax — 2m^/^keV. 

Strong differences (more than 10 per cent for luminosity, 
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5 per cent for the maximal energy) in the spectral energy 
distributions appear above m ~ 25 and rapidly increase with 
the mass accretion rate. This is considerably lower than the 
rhcT- ~ 100 Eddington limit. 

Here we assume that all the energy is released in the 
disc and that all the radiation is thermalized. First assump¬ 
tion is not strictly true because viscous stresses should be 
present eve rywhere throughout t he transonic fl ow (s ee sim¬ 
ulations bv IShafee et al.l feoOSl l: iPenna et al.l ll2012h l. The 
energy released in the inner parts of the flow contributes to 
the temperature of the gas and ultimately to comptonization 
of the soft radiation of the disc. We do not expect the op¬ 
tical depth of th e flow to true absorpt i on to be significantly 
large. Following iRvbicki fc LightmanI lll986ll . one can write 
down the free-free absorption cross-section for fully ionized 
hydrogen plasma as 

4e« 1 - 

where is electron density that may be substituted as 
M <? R m 

He — - —- (40) 

A-KHRu^nip GMgt H ^ ’ 

that implies the effective optical depth of 


^eff — '^z 



~ 5 X 10"® 


^2(yT')3/2 



(41) 


For a broad range of parameter values (E ^ T ^ Tmax) 
including the critical mass accretion rates estimated in this 
work, the effective optical depth is Tss 1. 


4 DISCUSSION AND CONCLUSIONS 

4.1 Limits of the standard disc approximation 


We considered the balance of forces acting upon a test parti¬ 
cle on the surface of a thick disc. Apart from radiation pres¬ 
sure due to electron scattering there are other forces that 
can induce launching winds from an accretion disc includ¬ 
ing magnet i c stre sses and radiation pressure in emission lines 
(see lPro^ ll2007l l for review about line-driven disc winds in 
active galactic nuclei) hence it is probable that in fact any 
disc is forming outflows even below the Eddington limit. 
However, the continuum (at least, electron-scattering) op¬ 
tical depth of the outflow may become considerably high 
only above the critical accret ion rate limit (re l evant optical 
depth estimates are given bv IPoutanen et aP l|2n07l l. equa¬ 
tions (29) and (31)). 

One of the main questions in the studies of supercritical 
accretion is the relative importance of energy advection and 
matter loss. The two limiting ca ses are the outflow-regulated 
accretion model introduced by IShakura fc Sunvaey Jl97^ 
and th e totally conservative slim disc model (see ISadowskil 
ll201lll for review) in which the larger part of thermal energy 
released by viscous stresses is advected inwards and finally 
lost in the black hole. Both models predict lower accretion 
efficiency (per unit mass in-flowing at large distance) and lu¬ 
minosity exceeding the Eddington limit by a logarithmic fac¬ 


tor. Similar logarithmic factor In ( | arises if one inte- 

\ Rin J 

grates the locally Eddington flux oc r~^ upon the surface of a 
disc with the outer and inner radii equal to Ro^t and Ri^. An¬ 
other effect usually neglected in accretion disc calculations 
but probably important for thick discs is radial radiation dif¬ 
fusion. This process tends to redistribute the r adiation flux 
according to von Zeipel rule dvon Zeipellll924ll F oc g (flux 
proportional to local gravity) that scales identically with the 
local Eddington limit. Hence there are different reasons for 
supercritical discs to attain the Teft oc asymptotic for 

the super-critical regime and the following logarithmic lu¬ 
minosity growth and “flat” spectrum {vLi, ~ const), and 
individual effects are difficult to disentangle. 

We find that transition to supercritical regime already 
appears for considerably large advection: only about 10 per 
cent of the energy released in the disc can escape as radia¬ 
tion and involved in a cceleration of dis c win ds. This justifies 
the us e of models like iLipunoval lll999ll and IPoutanen et al.l 
ll2007ll that include both outflows and advection. In such 
models, the mass accretion rate upon the black hole is con¬ 
siderably lower but still of the same order with the inflow 
rate at the outer rim of the disc. The overall radiative effi¬ 
ciency of mass accretion is probably close to 0.01—0.02 with 
respect to the mass absorbed by the black hole. 


4.2 Comparison to the slim disc model 

Our consideration still does not account for some other ef¬ 
fects that become important close to the limit. Firstly, the 
transonic nature of the thick disc solution that is incorpo¬ 
rated in the slim disc solution was not considered here. On 
the other hand, the slim-disc model does not account for 
the vertical structure but rather operates w ith vertically- 
averaged equations. In ISadowski et al] ll201lll . the vertical 
structure of a slim-disc was considered. However, the verti¬ 
cal structure itself was in this paper assumed to be decoupled 
from the radial structure of the disc, the vertical epicyclic 
frequency (corresponding to Qz or of the present work) 
was assumed to be constant with the vertical coordinate z, 
as well as all the velocity components, that means that the 
disc was in fact treated as thin. This allows to view our 
model as complementary to the slim disc. 

Another difference with the slim-disc approach is our 
unconventional treatment of advection effects. If advection is 
the primary term that breaks the local equilibrium between 
dissipated and radiated energy, and the advection effect is 
small, it is natural to describe advection as a delay term In 
the energy conservation equation: 

Q+-Q- = ivV)U, (42) 

where U is the local surface energy density, and and 
Q~ are the flux generated inside the disc and lost by its 
surface. One can approximate Q~ ~ U /tdiff and arrive to 
the approximate relation 



(43) 

or, in a linear approximation. 


Q^{r + Sr) ~ Q~{r), 

(44) 


where Sr is the radial correction introduced in section [Ml 
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Besides, this approach has the clear physical meaning that 
the photons generated at some radius drift downstream be¬ 
fore being radiated. Of course, this is a linear approach. 
For strong advection, the classical treatment of advection as 
storage of some fixed part of the locally generated heat may 
prove to be more relevant. In this sense, we also find our 
model complementary to the slim-disc model. 


4.3 Observational implications 

It is often adopted for simplicity that the properties of accre¬ 
tion discs change abruptly above the critical mass accretion 
rate. However, it seems that for a large span of dimension¬ 
less mass accretion rates, approximately between 1 and 100, 
accretion may be considered close to standard but already 
altered by thickness and advection. Besides, for rh > 1, the 
optically thick part of the disc extends further inwards that 
is extremely important in continuum fitting. Detailed con- 
tinuum fitting i s important for estimating bla ck hole spins 
llLi et al.l 1200^ : iNaravan fc McClintockI 120121 '). Black hole 
spins are measured in the high, or soft, stat es associated 
with h igher mass accretion rates with m ~ 1 (iMever et al.l 
[200 o 3). For m = 1, corrections for vertical gravity and ad¬ 
vection are unimportant, but rapidly increase with rh. In 
future, more comprehensive models they should be taken 
into account together with the effects of radiation transfer 
and scattering in the inner parts of the flow. Since the inner 
disc edge is effectively shifted closer to the event horizon one 
should expect ignoring the radial advection to overestimate 
the spin of the black hole. 

Bright quasars and QSO (quasi-stellar objects) could 
be another group of objects that exist close to the critical¬ 
ity limit. The puzzling five per cent accretion efficiency re- 
quir ed by some observati onal data such as quasar microlens- 
ing llMorgan et al.l[201(]ll may be connected to mild advec¬ 
tion effects lowering the radiative efficiency of QSO accretion 
discs. 


4.4 Conclusions 

We summarize that, due to stronger vertical gravity, transi¬ 
tion to super-Eddington accretion should appear at higher 
mass accretion rates. Close to the local Eddington limit, disc 
thickness becomes large and several effects neglected in the 
standard disc model become important. Here, we estimate 
the effects of vertical gravity dependence on vertical coor¬ 
dinate and the effects of non-local emission of the energy 
released in the disc. The latter is primarily connected to ra¬ 
dial advection. We find that including both effects simulta¬ 
neously tends to further increase the local Eddington limit, 
but the efficiency of accretion decreases to about one-two 
per cent for the critical mass accretion rate. More compre¬ 
hensive studies taking into account the non-one-dimensional 
character of thick accretion disc structure are needed to de¬ 
cide upon the relative importance of different effects such 
as advection and outflow formation in accretion discs ex¬ 
ceeding the Eddington limit. It may turn in fact that there 
is no unique Eddington limit, but accretion disc properties 
change much more gradually with mass accretion rate. 
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APPENDIX A: KEPLERIAN VON ZEIPEL 
CYLINDERS IN KERR METRIC 


A1 Basic equations 


Let us assume that dimensionless angular momentum I — 
—Utpjut and angular velocity u = jv!' depend on each 
other unambiguously, I = l{ui). Note that I here is the an¬ 
gular momentum per unit energy rather than per unit rest 
mass. This is the case for barotropic flows where gradients of 
pressure and density are collinear everywhere. This sentence 
is known as Poincare-Wav re theorem an d is proven both 
for classical mechanics fseelTassoull l|200(]lll and for General 
Relativity (iFriedman fc StergioulasI l2013ll. Sometimes it is 
to as von Zeipel theorem jAbramowic 3|1971I: 

199dl . Surfaces of constant I and w coincide 
1990h and are known as von Zeipel cylinders, 
however rely upon the given relation between 

Let us assume that the accretion disc is strictly Keple- 
rian in the equatorial plane: 


also referred 


[Chakrab^^ 


I Chakrabart 


Their shapes 
I and uj. 


U = UJk 


1 

-I- a ’ 


(Al) 


_ r^ — 2ay/r + e? 
^ 3/2 _ 2s/r a 


(A2) 


Expression (EU allows to exclude the radius r = 
and evaluate angular momentum as a function 

of LJ-. 


1(.) = + (A3) 

a;i/3 h -2(l-aw)^/®a;2/3j 

Angular velocity may be then found as a function of the two 
spatial coordinates by solving the equation: 


u. 

L{UI) — -;-, 

gtt + i^g^t 


(A4) 


where the left-hand side applies (lA^ . while the right-hand 
side is a consequence of velocity conversion relations = 
giptu^ + gipipu'^, ut = gttu^ + g^tu'^■ Relevant metric coeffi¬ 
cients: 


, 2 r 
(?tt — —1 H—2 ’ 


(A5) 


where: 


2 ar(l — p^) 

9vt ^2 

(A 6 ) 

dvv ~ p2 ^ ’ 

(A7) 

2 2 n 

fl = COS 6, 

(A 8 ) 

2 2,22 
p = r + a , 

(A9) 

^r^ + a^f-a^A{l-p^), 

(AlO) 

A = r^ + — 2r. 

(All) 


Finally, the main equation for angular velocity becomes 
(1 — (1 — 3aai) -I- — 2ar 


T/3 (^l-2(l-aa;)^/®a;2/3j 


- 2r 
1-p? 


(A12) 


-f 2aruj 


A2 Schwarzschild case 

If one sets a = 0, it is possible to obtain an analytical expres¬ 
sion for angular velocity as a function of r and p. Equation 
(IA12II becomes 

This equation is cubic both in r (useful for recovering the 
shapes of the cylinders) and in / = . 

/^-i/^+X = 0, (A14) 

r -2 

where X = — -;r—^. This equation is easily transformed 

2(1 — p-‘)r-^ 

into a depressed cubic equation R+pt+q = 0 by a substitute 
t = f- 1/6: 

In this equation, p — —1/12 < 0. Equation discriminant 
D = 4p^ -I- 27q^ < 0 for any r > 6 and p < v^, hence: 

t = i cos ^^acos (1 — 108X) — , (A16) 

/ 1X3/2 

and u) = t H— . Keplerian rotation field in the 

V 6/ 

Schwarzschild case is shown in Fig. (EH). 


A3 General Kerr case: solution and results 

Right-hand side of equation (IA12II is singular at small lj, but 
when searching for the solution for given r and p, one should 
be conscious about the singularity in the left-hand side of 
the equation, setting the maximal possible frequency: 

1 - 2 ( 1 - a;S/ax = 0, (A17) 

l^max-‘i^max + ^ (^18) 

a Sa 
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Figure Al. Keplerian von Zeipel cylinders for a Schwarzschild 
black hole. Red thick line is the applicability limit of the solution 
given by iIXTeli . 


Substituting cumax = t H-, one arrives again to a cubic 

3a 

equation: 


A- 


1 2 
27a3 




For any a 
4 




27a6 
expressed as 


l+(l-^a^ 
' 16 


0 , p 
2 


1 

3a2 


-TTT < 0: 


(A19) 
D = 


> 0 hence the solution may be 


= -T^-sign — 


3a 


8a 27a3 


X cosh ( iacosh ( — 
16 


a I 1 — 


—I 

27a2 J 


(A20) 


For 0 < w < tUmax, equation IIA12|l was solved numerically 
for tu for different r outside the last stable orbit radius and 

for |/r| < —r=. Isolines of resulting angular frequency for 

v 2 _ 

different Kerr parameter values are shown in Fig. IA2I 
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Figure A2. Keplerian von Zeipel cylinders for different Kerr parameters. Dotted curves mark the outer limit of the ergosphere. Cross 
is the position of the last stable orbit in the equatorial plane. 
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